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We present an experimental study of the infrared conductivity, transmission, and reflection of a 
gated bilayer graphene and their theoretical analysis within the Slonczewski-Weiss-McClure (SWMc) 
model. The infrared response is shown to be governed by the interplay of the interband and the 
intraband transitions among the four bands of the bilayer. The position of the main conductivity 
peak at the charge neutrality point is determined by the interlayer tunneling frequency. The shift 
of this peak as a function of the gate voltage gives information about less known parameters of 
the SWMc model, in particular, those responsible for the electron-hole and sublattice asymmetries. 
These parameter values are shown to be consistent with recent electronic structure calculations for 
the bilayer graphene and the SWMc parameters commonly used for the bulk graphite. 
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I. INTRODUCTION 

Since the monolayer graphene was isolated 1 - and shown 
to exhibit the quantum Hall effect 2 ^ a few years ago, ul- 
trathin carbon systems have attracted tremendous atten- 
tion^ Their electron properties are quite unique. Mono- 
layer graphene has a vanishing Fermi point at the Bril- 
louin zone corner and low energy quasiparticles with a 
linear spectrum, e(k) = ±i>|k|, which obey a massless 
Dirac equation. Here k is the deviation of the crystal 
momentum from the Brillouin zone corner (K point), 
v = (3/2)jQa/H is the quasiparticle velocity, 70 is the 
nearest- neighbor hopping parameter, and a — 1.42 A is 
the carbon-carbon distance. Graphene is the basic build- 
ing block of other types of carbon materials. Indeed, 
the first calculation of its band structure by Wallace 5 
was motivated by his studies of graphite. Extending 
that work, Slonczewski and Weiss,- McClurepii 7 - and oth- 
ers^ have developed the now commonly used Slonczewski- 
Weiss-McClure (SWMc) model for the low-energy elec- 
tron properties of graphite. This model, which is equiv- 
alent to a tight-binding model with seven parameters^ 
has proven to be a very useful analytical tool. It per- 
mitted theoretical calculations of a vast number of prop- 
erties of graphite, including its diamagnetic susceptibil- 
ity, de Haas- van Alfven effect, magnetooptical response, 
cyclotron resonance, and so on. These properties were 
actively studied experimentally until the late 70's and 
lead to accurate estimates of the principal SWMc pa- 
rameters, 70 through 73. Still, it proved challenging to 
unambigously determine the remaining three SWMc con- 
stants 74, 75, and A, which are measured in tens of meV. 
For illustration, in Table U we list inequivale nt p arame- 
ter sets from the latest original sources, Refs. Ildl and [TTI. 
Subsequently, the issue was further confounded by nu- 
merous misprints in reference books and reviews^. The 
density-functional theory c alculat ions ) 1 3 ' 14 ' 1 5 which nor- 



mally have accuracy of ~ 0.1 eV for quasiparticle disper- 
sion, have not yet settled this discrepancy. 

In view of the reinvigorated interest to graphene, it 
has become an important question to obtain the SWMc 
constants for a few layer graphene and also to compare 
them with those for bulk graphite. Thus, some differ- 
ence between the graphite and a graphene bilayer was 
recently reported, based on the analysis of Raman scat- 
tering.— Several ab initio calculations of these parame- 
ters for the bilayer have also been done i 17 ' 18 ' 19 ' 20 ' 21 Un- 
fortunately, they have not explicitly discussed the less 
accurately known SWMc parameters. 

The bilayer is a system intermediate between graphene 
and bulk graphite. Its lattice structure (for the case of 
the Bernal or AB stacking) is illustrated in Fig. [lja) . The 
corresponding band structure ; 23 ! 24 ' 25 shown in Fig.QJb), 
consists of four bands. These bands arise from splitting 
and hybridization of the Dirac cones of the individual lay- 
ers by the interlayer hopping matrix element 71 and by 
the electrostatic potential difference V between the two 
layers ! 23 ' 26 The latter can be controlled experimentally 
by varying the voltage V g of a nearby metallic gat o 27 ' 28 
or by doping^ This degree of tunability makes the bi- 
layer graphene an extremely interesting material for both 
fundamental study and applications. 

In this paper we show that 71, V4 = 74/70, and A 
can be directly extracted from the dynamical condictiv- 
ity measured in zero magnetic field. This is in contrast 
to the bulk graphite where determination of the SWMc 
constants was never straightforward and almost invari- 
ably required the use of strong magnetic fields. 

The dynamical conductivity er(f2) is determined by 
the six possible transitions among the four bands, see 
Fig- Die). They have energies of the order of a few 
1CT 1 eV, which is in the infrared optical range. Recently, 
experimental measurements of the infrared response of 
the bilayers have been carried out by our— and othe r 31 ' 32 
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TABLE I: The SWMc parameters (in eV) according to previous and present work. The numbers in parentheses are the reported 
accuracy of the trailing decimals. The "Exp" and "DFT" stand for experiment and density functional theory, respectively. 
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a L. M. Malard et al., Phys. Rev. B 76, 201401 (2007)i& 

b U. Min, B. R. Sahu, S. K. Banerjee, and A. H. MacDonald, Phys. Rev. B 75, 155115 (2007). 20 

C M. S. Dresselhaus and G. Dresselhaus, Adv. Phys. 30, 139 (1981)i£ 

d R. O. Dillon, I. L. Spain, and J. W. McClure, J. Phys. Chem. Solids 38, 635 (1977)ii 

e R. C. Tatar and S. Rabii, Phys. Rev. B 25, 4126 (1982)i2 

fj.-C. Charlier, X. Gonze, and J.-P. Michenaud, Phys. Rev. B 43, 4579 (1991)ii 
9M. Orlita et al., Phys. Rev. Lett. 100, 136403 (2008)^2, 
h A. Griineis et al., Phys. Rev. Lett. 100, 037601 (2008)i£ 

'This value cannot be very accurately found from our analysis and is instead adopted from the literature. 
^Physically irrelevant in the bilayer but should be set to zero for calculating A' from A. 
*Our estimate based on digitizing band dispersion graphs published in Refs. 18 19 20 21. 
'Absolute value only. 

""The negative sign (omitted in Ref. Il5ll is required for consistency with the conventional definition! of A. 



groups. Below we identify and explain the key findings 
of these experiments based on how different combina- 
tions of the interband transitions are either activated or 
suppressed by the Pauli exclusion principle. Our theory 
enables us to reach a quantitative agreement with the 
experiment using SWMc 70, 71, 74, and A, and also the 
phenomenological broadening constant T as adjustable 
parameters. The values of the SWMc parameters that 
give the best fit are given in the second column of Ta- 
ble |TJ Note that the next-nearest layer hopping parame- 
ters 72 and 75 are irrelevant for the bilayer. The parame- 
ter 73 cannot be reliably estimated from these particular 
experiments because it has an effect similar to the sim- 
ple broadening (r) in the range of carrier concentrations 
suitable for our analysis. 

Previous theoretical studies of the optical conductivity 
of bilayer graphen o 24 i 25 ' 33 i 34 ' 35 used a simplified model 
in which only 7 and 71 were taken into account. This 
model successfully explains the major features of cr(f2) 
as well as its dependence on the gate voltage V g , and we 
qualitatively summarize it as follows. Conduction and 
valence bands are symmetric. In the absence of the elec- 
trostatic potential difference V between the layers the 
two conduction (valence) bands have the same shape and 
are shifted by 71. Except the range of very small mo- 
menta k, their shape remain nearly identical even in the 
presence of a finite V. As a result, there is a high op- 
tical density of states for transitions between the two 
pairs of bands at frequency 71 /h, which gives rise to a 
sharp peak in the real part of the conductivity Recr(f2) 
at fl = 71/fi w 3200cm" 1 (using 71 = 0.40 eV). Other 



transitions give more gradually varying contributions to 
Recr(f2), eventually leading to the asymptotic "univer- 
sal" vsAue&2&2M2£^2&ZL a = e 2 /2h at high frequency 
(which is twice the value for the monolayer—). Finally, 
in real graphene systems the conductivity features are 
never sharp because of a finite lifetime due to, e.g., dis- 
order scattering. This broadens the peaks and can also 
merge together several features that are close in energy, 
see Fig. H 

Our recent infrared experiments 3 ^ as well as measure- 
ments by another group^ 2 . have largely confirmed this 
picture but also found features that cannot be explained 
within this simple model. In particular, the conductivity 
peaks on the electron and the hole sides are displaced 
in energy from 71 by about 10% in opposite directions. 
[Electron and hole doping is identified with, respectively, 
positive and negative SV = V g — Vcn, where Vcn is the 
gate voltage at which the bilayer is tuned to the charge- 
neutrality (CN) point.] 

In order to investigate the origin of these features 
in this paper we carry out a combined experimental- 
theoretical study of the infrared response of a bilayer 
graphene. We attribute the observed electron-hole asym- 
metry to the effect of 74 and A. We find that including 
these parameters is essential for a more accurate discus- 
sion of <r(f2) of the bilayer. Besides differences in the 
optical response, 74 and A also make effective masses 
for electrons and holes unequal* 3 ^ in agreement with the 
findings from the Raman scattering^ 

In our experiments, we have measured the optical re- 
flection R(ft, V g ) and transmission T(f2, V g ) as a function 
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of the frequency fl and the gate voltage V g . Prom R and 
T we extracted the real and imaginary part of the con- 
ductivity using a commercial software package. Some of 
these experimental results were reported previously^ 

In this paper we present more extensive experimental 
data and we also compute the same three quantities - 
a, i?, and T — theoretically. The calculation requires 
accounting for the interplay of several physical phenom- 
ena: (a) electrostatic charging of the layers (b) their dy- 
namical conductivity, (c) disorder, and (d) the optical 
properties of the environment (sample, substrate, and 
the gate). Each of these ingredients has been studied in 
the pastM£&2LM£M£d2dL Here we carry out all these 
calculations in a single paper albeit we include disorder 
broadening in a very simple way. This enables us to di- 
rectly compare our theoretical results with the measure- 
ments. 

The remainder of the paper is organized as follows. In 
Sec.|TT]we summarize our results. Theoretical derivation 
is outlined in Sec. IIIII Section HVl contains comparison of 
the theory and experiment, discussion, and conclusions. 
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FIG. 1: (Color online) (a) Crystal structure of the graphene 
bilayer with the relevant SWMc hopping parameters shown, 
(b) Band structure of a biased bilayer (lines), which can be 
considered as hybridization of two shifted Dirac cones (dots) . 
Numbers on the right label the four bands, (c) Examples 
of the allowed optical transitions for the chemical potential 
indicated by the dashed line. Occupied states are shown by 
the thicker lines. The dots and the arrows mark the initial and 
the final states, respectively, of the transitions that produce 
features at frequencies Ej, j = 1, 2, . . . , 6 in Fig. [2ja) below. 
Eq is the intraband transition (Drude peak). 
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FIG. 2: (a) Real and (b) imaginary part of conductivity in 
units of (To = e 2 /h for the gate voltage SV = — 100 V. The 
solid curves are for broadening V = O.O271. The dashed curve 
is for r = O.OO271. 



Some calculational details are relegated to the Appendix. 



II. RESULTS 

To measure the optical response of the bilayer we em- 
ployed synchrotron infrared radiation, as described pre- 
viousry. 30 ' 42 Understandably, the two-atom thick sample 
has a rather small optical signal. The quantity which 
can be extracted most reliably from the current exper- 
iments is the relative transmission T(f2, V g )/T(Q, Vcn) 
and reflection i?(fi, V g )/R(tl, Vcn)- All measurements 
were done at the temperature of 45 K. The data for the 
largest |<$V| = \V g — VqnI are depicted in Fig. [3) The main 
feature in the relative transmission spectra is a small but 
clearly visible dip around f2 = 3200 cm -1 . Away from 
the dip, the relative transmission is slightly higher than 
unity. The relative reflection spectra are characterized by 
a dip-peak structure. Transmission and reflection spectra 
are asymmetic between positive and negative 8V, which 
correspond, respectively, to doping of electrons and holes 
in bilayer graphene. 

From the transmission and reflection data, we ex- 
tracted the optical conductivity . 30 ! 42 ' 43 The dominant 
feature in the conductivity spectra is a strong peak at 
il w 3200 cnT' 1 , see Fig. Ht^c) . Below the main peak, we 
observed a broadened threshold feature, which shifts sys- 
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tematically with SV. The most intriguing observation is 
again the electron-hole asymmetry in the optical conduc- 
tivity. For instance, the frequencies of the main peak in 
Recr(f2) and its voltage dependence are noticeably dif- 
ferent for electrons and holes, see Fig. [He). Also, while 
the peak is quite symmetric at large positive voltages, at 
high negative SV , it is not. The most probable reason is 
the existence of a secondary peak at a slightly larger fl, 
see below. 

On the theory side, we calculated c, T, and R, using 
the SWMc constants and T as adjustable parameters. 
Results for the conductivity are shown in Fig. E^b). The 
reflection and transmission are plotted in Fig. [3l The 
calculational parameters were adjusted to reproduce the 
frequency positions and widths of the main features of the 
experimental data. Interestingly, in this way of fitting, 
it was not possible to achieve an equally good agreement 
for the vertical scale of the observed features. Still their 
qualitative trend as a function of SV is reproduced well. 

Both in experiment and in calculations the carrier con- 
centrations are always smaller than the characteristic 
value tiq given by 



n 



h 2 v 2 



3.7 x 10 13 cm" 2 . 



(1) 



Here and below we assume that 70 = 3.0 eV, which cor- 
responds to v = (3/2)7oa/7i = 1.0 x 10 8 cm/s. (Based 
on other results in the literature, this value should be 
accurate to about 10%.) At concentrations \n\ < no 
the high energy bands 1 and 4 have no free carriers and 
Rea(fi) has a pronounced peak at fi w 3200cm -1 . As 
explained above, this feature corresponds to transition 
between band pairs that are nearly parallel: bands 3 and 
4 for (i > or bands 1 and 2 for /x < 0, see Fig. [1] 

The evolution of the infrared response with V g can be 
understood as follows. As the gate voltage deviates fur- 
ther away from Vcnj the electron concentration 



n = C b SV/e 



and the chemical potential /i increase by the absolute 
value. Here C b is the capacitance between the bilayer 
and the gate. As a result of an increased |n|, the peak be- 
come more pronounced. Simultaneously, near the higher 
frequency side of the peak a depletion of conductivity 
develops. One can say that the optical weight is increas- 
ingly transfered from the high frequencies to the 71 peak. 
Larger conductivity is directly associated with decreased 
optical transmission. Therefore one observes an increas- 
ing dip in the transmission near 71 and a higher trans- 
mission at higher f2, see Fig. [3l Similar features appear 
in the reflection but they are more difficult to interpret 
as they are also affected by lm<r(Ct). 

Very important for our analysis are the aforementioned 
small shifts in the position of the 71 peak as a function of 
SV. Within the SWMc model, their origin is as follows. 
In the absence of broadening, the peak arises from the 
absorption in the range of frequencies, Ei < Kl < £3, see 
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FIG. 3: (Color online) Relative transmission: (a) theory (b) 
experiment. Relative reflection: (c) theory (d) experiment. 
The solid line is for electrons, 8V « +80 V. The dashed line 
is for holes, SV ~ —80 V. The experimental uncertainties 
are ~ 0.002 (0.2%) at Q near 3000 cm -1 and ~ 0.5% at high 
frequency. 



Figs. Q] and [2j Since the optical weight at £3 is higher, the 
conductivity peak occurs at energy £"3. However, if the 
broadening is large enough, the optical weight becomes 
distributed more uniformly, and the peak position moves 
to the midpoint of E2 and £3, see Fig. [5] Energies £2 
and £3 themselves vary with the gate voltage (or n) . For 
positive SV (positive n) , E 2 = E% is the energy difference 



(2) between the bands 3 and 4 at k = 0. The energy £3 



Et is the corresponding difference at k = hp, where 



kp = sign (n) \/7r|n| 



(3) 



is the Fermi momentum. For SV < we denote Ei and 
E3 by, respectively, E^ and E^ and they are computed 
using the bands 1 and 2 instead of 3 and 4. 

From the band structure ; 24 ' 25 we can find the following 
approximate expressions valid for n <C no: 



(4) 




± A =F 2 (2«47i 



(5) 



Here V = V(n) as well as the chemical potential u = 
n(n) are determined self-consistently by the electrostatics 
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FIG. 4: (Color online) (a), (b) Theoretical and (c) experimental results for the conductivity Re<r, in units of cto = 4e 2 //i, as 
a function of frequency Q. The deviation SV of the gate voltage from the charge neutrality point is indicated next to each 
curve. For clarity, the curves are offset vertically by 0.5<7o from one another. The SWMc parameters for plot (b) are given in 
Table U In (a) they are the same, except 73 is set to zero. The dashed curves superimposed on the SV = +50 V (— 50 V) traces 
in (a) and (b) are the arithmetic means of all the positive (negative) SV curves. Their significance is discussed in Sec. IIVI The 
estimated uncertainty of the measured Recr is 0.125<to at f2 ~ 8000 cm -1 and 0.0625<jo at f2 ~ 3000 cm - . 



of the system^ see Sec. IIIII These equations indicate 
that the parameters primarily responsible for electron- 
hole asymmetry are 74 and A. 

Parameter A is the difference of the on-site electron 
energies of the A and the B sites&2 [the stacked and un- 
stacked sublattices, respectively, see Fig. Ufa)]. It has 
two effects: first, it lifts the k = energy for bands 1 
and 4; second, it adds a k dependent perturbation to the 
two band dispersion. Parameter V4 = W74/70 of dimen- 
sion of velocity characterizes hopping between a stacked 
atom and its the three unstacked neighbors of its stack- 
ing partner. It also introduces difference between the 
valence and conduction bands. To the leading order in k, 
this hopping shifts the two middle bands (2 and 3) up- 
ward by a term proportional to v^k and shifts the two 
outer bands (1 and 4) downward by the same amount. 
These effects of A and V4 are illustrated in Fig. [6l 

Additional electron-hole asymmetry can in principle 
come from extrinsic sources, e.g., charged impurities that 
can be present on or between the layers. Besides creat- 
ing a finite Vqn, these charges also move V — point 
away from the charge neutrality point n = 0. To the first 
approximation,— this introduces an offset of the inter- 
layer bias: V(n) — ► V(n) + Vo- However, our calculations 
suggest that for reasonable Vq this effect has a smaller 



influence on the electron-hole asymmetry of the optical 
response than A and 74. 

Based on the above discussion, we can predict quali- 
tatively how the position of the main conductivity peak 
should vary as a function of SV. For example, on the 
electron side, and for V4 > 0, the peak should move to 
lower frequencies as SV increases. Alternatively, this can 
be seen from Fig. [6] the top two bands move closer to 
each other as k = kp increases. 

For the quantitative analysis, we use a full numerical 
calculation of a and T, which is discussed in Sec. IIIII bc- 
low. It demonstrates that for the case of small T the en- 
ergy E3 is indeed in a good agreement with the computed 
peak position Qq. However, the broadening observed in 
experiment o 30 ^ 31 ' 32 is appreciable, in which case the for- 
mula Slo = (£2 + Es)/2 is more accurate. Of course, 
for fairly large T other nearby transitions, £1 and £4, 
start to influence the lineshape of the main peak. This 
is especially noticeable on the hole side, where the £4- 
peak is right next to the main one. In the calculations 
this two-peak structure is unmistakable, see Fig. Sfb). 
In the experiment, where the main peak is for some rea- 
son strongly enhanced compared to the calculation, the 
£4 peak is somewhat disguised. As pointed out by Kuz- 
menkOf2£ the difference between £4 and £2 can in prin- 
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FIG. 5: Position of the 71 peak vs. gate voltage for the two 
values of the broadening: (a) T = O.O271 and (b) F = O.OO271. 
The solid lines are our numerical results from the conductiv- 
ity; the thick lines are from the relative transmission. The 
dashed lines show E3 and (E2 + Ez)/2 in the cases (a) and 
(b), respectively. The SWMc parameters used in the calcula- 
tion are listed in the first column of Table [T] except in (a) 73 
is set to zero. The symbols are the peak positions determined 
from the measured conductivity (squares) and transmission 
(triangles) . 



ciple provide a direct spectroscopic measurement of the 
energy gap V. 

For detailed comparison with experiment we use our 
numerical results rather than Eqs. @ and (JSJ). Fitting 
them to the data, see Fig. [5j we have obtained estimates 
of 71, 74, and A listed in Table Q] This fitting procedure 
proved to be very straghtforward. For example, A is de- 
termined mostly by the splitting of the peak positions 
on the electron and the holes sides of the charge neutral- 
ity point. Parameter 71 is essentially the average of the 
two. Finally, 74 controls the slope of the Qo(V g ) curves 
away from Vcn- Therefore, all these parameters can be 
uniquely determined. 

In Table U we also list SWMc values suggested in prior 
literature. They mainly agree with ours for the principal 
SWMc parameters 70 and 71 but show some deviations 
for the more subtle quantities 74 and A we have been 
discussing here. Possible reasons for these differences are 
given in Sec. IIVI 




FIG. 6: The effect of 74 and A on the band structure. Pa- 
rameter A raises the bands 1 and 4. The interlayer neighbor 
hopping term 74 gives a contribution quadratic in k opposite 
in sign for the conduction and the valence bands. The solid 
(dashed) lines are the bands with positive (zero) value of 74. 

III. DERIVATION 
A. Band structure 

The bilayer is two monolayers stacked together, see 
Fig. [Ha). I n the bulk graphite the preferential stacking 
is the AB (Bernal) one, such that only one sublattice of 
each layer is bonded to each other. In order to achieve 
agreement with experiments,— we have to assume that 
in the bilayer the stacking is the same. We use the ba- 
sis {"Fai, ^bi, ^B2, ^42}, where the letter stands for the 
sublattice label and the number represents the layer in- 
dex. In this basis the SWMc tight-binding Hamiltonian 
for the bilayer becomes^ 5 - 

/-T + A 7l -V>*\ 

TT V ~\ -V 4 (j)* V 3 (j) , , 

71 -V49 "2 + A 

\ -V 4 (f) V 3 (f>* \ J 

where = —i{k x + ik y ) and (k x ,k y ) is the deviation of 
the quasimomentum from the K point. 

Given V, it is easy to obtain the four band energies 
e a (k) and the corresponding eigenstates |a,k) numeri- 
cally. However, as mentioned in Sec. [TTJ, V should be 
determined self-consistently as a function of V g , or equiv- 
alently, the total carrier concentration n. The algorithm 
for doing so is given next. 

B. Electrostatics 

As discussed in the literature ; 2 ^ 40 the electric field of 
the gate has two major effects on the bilayer graphene. 
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FIG. 7: (a) Interlayer bias V as a function of total density 
n. Three sets of curves correspond to (from top to bottom) 
Vo = O.I71, 0, and — O.I71. The dashed lines are computed 
from Eq. (I12[) . (b) Chemical potential vs. n for Vo = — O.I71. 



First, it modifies the bands by introducing a potential 
difference between the layers and as a consequence opens 
up the energy gap. Second, it induces charge carriers. 
Electric field of the charged impurities can play a similar 
role: it creates a layer asymmetry Vo and opens a gap at 
the charge neutral point much like an external gate. But 
the more important effect of the impurities is presumably 
the broadening of the electron energy states, which we de- 
scribe by a phenomenological constant T. For example, 
if the impurities are distributed symmetrically between 
the two layers, then Vo is zero but T is still finite. We as- 
sume r to be real and independent of energy, momentum, 
or a band index. This is certainly a very rudimentary 
treatment of disorder compared to, e.g., self-consistent 
schemes i 25 ' 36 i 44 However, since the source of disorder in 
graphene is still debated, we think that this simple ap- 
proach is adequate for our purposes as long as T is treated 
as another adjustable parameter. 

To compute V(n) and fi(n) we set upa system of equa- 
tions similar to those in Refs.l4fjlandl24l. These equations 
capture the dominant Hartree term of the interaction 
but neglect exchange and correlation energies^ The first 



equation is [cf. Eq. ([2])] 

n = n t + rib = CbSV/e , 



(7) 



where nt and n& are the carrier concentrations of the top 
and bottom layers, and Cb is the capacitance to the gate. 
Second, the electrostatic potential difference between the 
two layers V is given by 



47re 2 

V = (n t - n b )c , 



(8) 



where k is the dielectric constant and Co is the distance 
between the layers. Next, the Hamiltonian and hence the 
wavefunction and the layer density n t and n& depend on 
V . Therefore the quantities V, n t , and n& must be solved 
for self-consistently. If the broadening T is neglected, this 
can be done analytically in the limit V, [i <C 71 , which 
gives V ~ V(n, Vo), where 



,,24,40 



V(n, Vo) = 



X-yx + V 



A- 



X 



I In |*l 



A = -, (9) 



no is defined by Eq. fl}, an d A = e 2 cono/(irKji) is the 
dimensionless strength of the interlayer screening. Using 
the typical parameter values, one estimates 4 ^ A ~ 1, and 
so the interlayer screening is significant) 20 ' 40 

For experimentally relevant broadening T ~ O.O271, 
the approximation leading to Eq. ^ is no longer accu- 
rate. Therefore, we computed the dependence of n* and 
rib on V numerically as follows. We first define the re- 
tarded Green's function G R by the analytic continuation 
G R (e) = G (e — > e + iT) of the following expression 



a=l ■ J 



(10) 



Then we compute rit from 



n t = 



i 1 



')]> (11) 



using numerical quadrature. Similarly, the formula for 
rib is obtained by replacing Gn + G22 with G33 + G44. 

The system of nonlinear equations (0, ©, and (fTTj) 
is solved by an iterative procedure. For a given chem- 
ical potential /x we start from some initial guess on V. 
Then we diagonalize the Hamiltonian and compute G^, 
n t , and rib. Substituting them into Eq. ([8]), we get the 
value of V for the next iteration. (Actually, we use not 
this value directly but a certain linear combination of the 
new and old V to achieve convergence.) The iterations 
terminate when the values of V changes by less than a 
desired relative accuracy (typically, 10~ 5 ). The results of 
these calculations are in a good agreement with Eq. (J9]) 
for r = 0, and so are not shown. On the other hand, 
the results for V = O.O271, which are plotted in Fig. [71 
appreciably deviate from Eq. ^ . The agreement greatly 
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improves (see Fig. [7]) if instead of Eq. ([9]) we use, on 
heuristic grounds, the following formula: 

V(n)=V(n*,V o )-V(n r ,0), (12) 

n* = sign(n)wn 2 + n% , np = sign(n) . (13) 

V 7T7l 

C. Dynamical conductivity 

The above procedure enables us to compute V and n 
for a given chemical potential /i. With the former deter- 



mining the Hamiltonian and therefore its eigenstates, and 
the latter determining their occupancy, we can now com- 
pute the dynamical conductivity by the Kubo formula^ 



n«(ft)-n«(o) 

1 1 + lO 



where the polarization operator U R x (Cl) is given by 



n* x (n)=ig ¥ J'J -Tr{v x [G R (k,e)-G A (k,e)]v x [G R (k 7 e + n) + G A (k,e-n)}} . (15) 

— oo 

In this equation g = 4 is the spin-valley degeneracy of graphene, v x = frr 1 dH/ dk x is the velocity operator, and G R - A 
at the retarded and the advanced Green's functions. Assuming again that the broadening is momentum and energy 
independent, these functions are obtained by the analytic continuation of G in Eq. (flU|) : G RA (e) = G (e — > e ± iT). 
After some algebra, we find 



II* (Q) = i 5 (|) 2 f-^£|M Q/3 (k)| 2 tK[sfi{k)-iTZ,e a {k)-(ir + Sl)C], 



(16) 



r 



where M a /3 (k) = (a, k jv^j /3, k) are the transition ma- 
trix elements and function K is defined by 



K{z\, z 2 ) 



ln(/x - zi) - ln(/i - z 2 ) 
2w(z 1 - z 2 ) 



(17) 



with the branch cut for lnz taken to be (— oo,0]. 

For vanishing V and T the conductivity can be com- 
puted in the closed form, see Appendix [B] For other 
cases, we evaluated it numerically. The results are shown 
in Figs. [5]and[U To demonstrate agreement with pre- 
vious theoretical calculations , 24 i 25 i 33 i 35 we present a(Q) 
computed for a very small broadening T in Fig. O In 
this case one can easily identify all six transitions. As 
explained above, the sharp features at fi « 3200 cm -1 
are due to the high optical density of states at energies 
E 2 < Kl < E3. The other prominent feature at SI = 
is the intraband Drude peak. (Its height is related to 
the transport mobility.) In Fig. [4] the calculation is done 
for much larger T to match the experimental data. This 
Figure has been discussed in detail in Sec. ITT1 



IV. DISCUSSION 

In this paper we presented a joint experimental and 
theoretical study of the infrared response of a bilayer 
graphene. Our results demonstrate a complex interplay 



among various interband transitions and their disorder- 
induced broadening. Nevertheless, by means of a care- 
ful analysis, we have been able to explain the majority 
of the observed features within the conventional SWMc 
model. The corresponding SWMc parameters are given 
in Table HI together with their estimated uncertainties. 
In particular, our 71 should have a very high accuracy: 
about 100cm -1 , i.e., 3%. The uncertainty in 71 comes 
predominantly from an unknown systematic error that 
we make by neglecting the renormalization of the spec- 
trum by scattering processes. Since we assume that the 
imaginary part T ss 65 cm -1 of the electron self-energy 
due to scattering is constant, its real part has to van- 
ish by the Kramers-Kronig relations. In fact, this real 
part, which is generally finite^ can shift the observed 
transition frequencies by an amount that scales with T. 

Let us now compare our SWMc parameters with those 
found in previous work on bilayers and bulk graphite. 
For the bilayer case there is at present only one other 
experimental determination^ of Jj's. From Table U we 
see that the difference between our and their values is 
primarily in 71. Actually, our SWMc parameters can 
describe the Raman data equally well^ as those given 
in Ref. [T5 . Our parameter values have smaller estimated 
errors and should be considered more accurate. 

In comparison with bulk graphite, the strongest dis- 
crepancy is in the value of A. The difference is signifi- 
cantly larger than the uncertainty of A grap hite quoted in 
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the earl y 10 ! 11 and the recent experimental work^ which 
makes a strong case that Abu aye r differs from Ag^^ite 
both in sign and in magnitude. To judge the true signif- 
icance of this result, one should recall that the physical 
meaning of Allayer is the difference in the onsite energies 
of the A and B sublattices.°& However, in graphite the role 
of the same quantity is played not by A but by the linear 
combination 9 

A graphitc = A graphitc - 72 + 75 • (18) 

For the sake of convenience, let us set 72 = 75 = in the 
bilayer, so that the A-B energy difference is equal to A' 
in both materials. Taking the most commonly usecU^ pa- 
rameter values for graphite, we arrive at the remarkable 
empirical relation 

A graphite « 37 meV « 2A bilaycr , (19) 

which is much easier to interpret. Indeed, the physical 
origin of A' is the short-range (exponentially decaying 
with distance) repulsion due to exchange and correlation 
effects between the electron states of the stacked atoms. 
(Neither Coulomb nor even the van der Waals interaction 
have short enough range to effectively discriminate be- 
tween the two sublatticeS ) 47 i 48 given the relatively large 
interlayer distance.) Since in the bilayer each A atom has 
a single stacking partner while in the Bernal graphite it 
has two of them, Eq. (|19p is exactly what one would 
expect. More precisely, it is expected if the interlayer 
distance in the bilayer and in the graphite are nearly the 
same. The validity of Eq. (|19p can be considered an ex- 
perimental evidence that this is indeed so. 

Another SWMc constant, which may seem to be dif- 
ferent in the bilayer and the bulk graphite is 74. As 
mentioned in Sec. U this is one of the parameters that in 
the past have been difficult to determine very accurately. 
Our estimate of 74 can be defended on the grounds that 
(i) it agrees with the Raman experiments 1 ^ and (ii) it 
is comparable to the accepted value of 73. These two 
parameters describe hopping between pairs of atoms at 
equal distances in the lattice, see Fig. QJa), and theoret- 
ically are not expected to be vastly different from each 
other. Large difference of 74 between the bilayer and the 
bulk graphite is not expected either. Indeed, even when 
they disagree about the order of magnitude (or sign) of 
A, all electronic structure calculations to date find that 
74 ~ 73 and are of the same order of magnitude in the 
two systems, see Table fl] 

Parameter 73 itself cannot be reliably extracted from 
the experimental data^ - we analyzed here. At the rele- 
vant carrier concentrations the main effect of 73 is to pro- 
duce a weak trigonal warping of the band dispersion fi - 
This warping averages out over the Fermi surface, and 
so has an effect similar to the broadening T: it makes 
the 71 conductivity peak more symmetric and shifts it 
towards the midpoint of E% and -E3, i.e., to slightly lower 
frequencies, cf. Figs. IUa) and (b). Thus, it is difficult 
to separate the effect of 73 from the broadening due to 
disorder. 



Regarding the latter, the dc mobility that we find 
from our numerically computed cr(0) using T — O.O271 ss 
8 meV is fl rj 3900 cm 2 / Vs. This is close to the trans- 
port mobility typical for bilayer graphene, supporting our 
interpretation that T arises mainly due to disorder. 

Concluding the paper, we wish to draw attention to 
several features of the experimental data that are not 
accounted for by our model. One of them is an unex- 
pectedly large amount of the optical weight in a range 
of frequencies below the 71 peak. It is present be- 
tween the Drude peak and 2/i, i.e., twice the chemi- 
cal potential. For the chosen T, our calculation pre- 
dicts Recr(f2) ~ 0.02e 2 /h at such f2, see Fig. [3J whereas 
the measured value is several times larger^ This extra 
weight is present also in the monolayer graphene, in the 
same range of frequencies.— A related issue is a very 
gradual rise of Re a(fl) around the point £1 = 2/i com- 
pared to a sharp threshold expected theoretically. These 
features can be in part due to electron-phonon interac- 
tioni 4 - or midgap state s 44 ! 49 but other effects seem to be 
involved as well. 

One very simple explanation would be to attribute 
both the broadening of the = 2/i threshold and the 
extra weight at Q < 2/i to long-range density inhomo- 
geneities in the sample. They can be caused by charge 
impurities and remnants of the photoresist used in the 
sample processing. The presence of such inhomogeneities 
would modulate the local chemical potential, and so in 
the infrared response one would see a certain average of 
the a (f2) taken at different SV. We illustrate this ar- 
gument by calculating the arithmetic mean of er(fi)'s 
for positive (negative) SV and superimposing the re- 
sults (shown by the dashed lines) on the a (f2) traces for 
SV = +50 V (-50 V) in Fig.gjb). Such averaged conduc- 
tivities indeed resemble the experimental data [Fig.fJJh)] 
more faithfully. 

Another discrepancy between the experiment and the 
present theory is the lineshape of the 71-peak. By vary- 
ing r, we can fit either the width or the height of the 
peak but not both. For example, in Fig. [4j where we 
chose to fit the width, the measured height of the peak is 
sometimes nearly twice larger than the theory predicts. 
The extra optical weight of the peaks appears to have 
been transferred from their high-frequency sides, which 
are suppressed in experiment compared to the calcula- 
tions. These lineshape differences are significant enough 
to make us think that some essential physics is still miss- 
ing in the simple single-particle picture presented in this 
paper. We speculate that including many-body effects 
may be truly necessary for bringing theory and experi- 
ment to better agreement. 
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APPENDIX A: REFLECTION AND 
TRANSMISSION 

To compute the transmission coefficient T and the re- 
flection coefficient R we follow the standard procedure^ 
In general, the result depends on the angle of incidence 
and on the polarization of light. Abergel and Fal'kc 51 
derived the formulas for R and T for the S'-polarization 
where the electric field is perpendicular to the plane of 
incidence (and parallel to the sample surface). We repro- 
duce them here with a slight change in notation: 



R 



T = 



C ni cos 6\ — D [cos 6a — Ana] 
C m cos Q\ + D [cos 9q + Ana] 

2 cos 9q n\ cos B\ Ti2 cos 82 
61+ D [cos 9 + Ana] 



C m cos 1 



where A, B, C, and D are given by 

sin 82 
> cos 82 



(Al) 



A = cos 62 sin 82 + iri2 cos 0q cos 82 , 
B = i cos #2 cos 82 + ri2 cos 9q sin 82 , 
C = Ani cos 9 2 sin 8± + iB n 2 cos 9\ cos <5i , 
D = iAni cos 62 cos 8\ + B n 2 cos 9\ sin 8\ . 



(A2) 



In Eqs. (jAip and (|A2j) . the index j = 0,1,2 represents 
vacuum, Si02, and Si layers respectively, rij are the index 
of refraction of each layer, and 9j are the angles the light 
ray makes with the surface normal in each layer. They 
satisfy Snell's law rij sui9j = const. Finally, 8j — kLj/rij 
is the phase the light picks up as it makes one pass across 
the layer of thickness Lj . 

For the other, P-polarization, where the electric field 
is not exactly parallel to the surface of the sample, we 
find a different expression: 



R = 
T = 



C n\ cos 9o — D cos 9\ (1 — Ana cos 0q) 
C ni cos 9o + D cos 9\ [1 + A-ko cos 8q] 

—2 cos 9q rii cos 9 1 TI2 cos 62 
C rii cos 9q + D cos 9i [1 + Ana cos 9q] 



(A3) 



in the infrared beam. Thus, we decided not to include it 
in the analysis and do all the calculation assuming the 
5-polarization only. 

APPENDIX B: CONDUCTIVITY OF AN 
UNBIASED BILAYER AT VANISHING 
BROADENING 

The conductivity for the case T = V = was com- 
puted previously in Refs. [H and [HI. In our attempt to 
reproduce their formula we discovered that it contains a 
typographical sign error^ For future reference, we give 
the corrected expression below. 

In the limit of zero broadening, r — > 0, Eqs. fT4" )) -(fl"7] ) 
reduce to the following expression for the conductivity: 



da; n\M, 



a/3 I 



(Bl) 

where P means principal value and the integration vari- 
able uj = \e a — Ea\ is the energy difference between two 
states. The sum in Eq. (IBlj) is over all values of momen- 
tum kj{uS) of which two states differing in energy u> exist. 
For V = where the the matrix elements M a p take a 
simple form, the integration over u> in Eq. (|B1|) can be 
done analytically. The result can be written as a sum of 
three terms: 



a{n) 

cto 



a (SI) + a yi (SI) + ct 27i (SI) , (B2) 



where 00 = e 2 /ft is the unit of conductivity, ao is contri- 
bution from transitions between bands 2 and 3 that turn 
on at SI = 0, <7 Tl is contribution from transitions between 
bands 1 and 3 and bands 2 and 4 that turn on at SI = 71 , 
(J2 7l is contribution from transition between bands 1 and 
4 that turn on at SI = 271. They are given by 
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7i 



02 7l 



1 SI + 2 7 i 

n 

10-27 



n 7 i 



2 Q(Sl 



7T 7 2 

i 

7T 



o 2 

271 
si 



In 



SI 

7i 

± 
SI 2 



(B3a) 



In 



71 



SI 



li -si 
-2 7 f 



(B3b) 



x In 



2 SI 

2 7 i +SI 



ie(o-2 7l )--(^ > 

71 tt\2 SI 2 - 7f 

4 7l 2 - SI 2 



2 7 i - SI 



1 Sin 
'2 SI 2 



■In 



7i 



7i 



For this polarization the conductivity enters R and T 
multiplied by the cosine of the angle of incidence, i.e., its 
effect is reduced. In our experiments, we typically have 
#0 ~ 30°, and so this reduction is quite small. Its role is 

further diminished by the presence of both polarizations where, for ease of notation, SI stands for hSl and g = 4. 



(B3c) 
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